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Abstract: The coastal zones with tidal mudflats have the most development potential for agriculture and industry in China. It is 
therefore vital to study the spatial patterns of soil heavy metals and to identify the intrinsic relationship between the spatial 
patterns of soil heavy metals and the natural and human driving factors in coastal regions. Using classical statistical and 
geostatistical methods, this study investigated the effect of land use types on the contents of total and bioavailable heavy 
metals (Pb, Cr, Cd and As) in surface soil not only to delineate the spatial distribution of these forms of soil heavy metals, but 
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also to determine the relationship between contents of bioavailable heavy metals and basic soil physical and chemical 
properties. The study was conducted in a coastal area in Jiangsu Province in China — a rapidly developing industrial and 
agricultural region. Results indicated that soil environmental quality across the study area was generally good. The contents of 
soil Pb, Cr, Cd and As were lower than the accepted environmental quality evaluation standard (the first class), although there 
was an increasing trend of accumulation of each heavy metal in recent years. Also environmental risk and ecological toxicity 
of soil bioavailable Pb, Cr, Cd and As in the study area was not significant. Land use type affected total contents and 
bioavailable contents of Pb, Cr and Cd to different degrees. Greenhouse facility soils and greenery park soils had the highest 
contents of total and bioavailable Pb. Also soil contents of total and bioavailable Cr were high in residential and industrial 
areas, whereas greenhouse facility soils had the lowest total and bioavailable Cd contents. Land use patterns had little impact 
on total and bioavailable contents of As. In addition to land use type, biological absorption was another important factor 
affecting the contents of soil heavy metals in the study area. Total soil Pb, Cr and Cd and bioavailable Pb and Cr had 
significant directional trend across the study area. Spatial pattern of the contents of total and bioavailable soil heavy metals 
was generally influenced by various factors at different scales. The whole spatial distribution of soil heavy metal was a 
strip-like pattern due to large-scale tidal deposits. Patched local distributions were mostly controlled by small-range factors 
such as industrial emissions and human activity. Soil bioavailable Pb, Cr and Cd exhibited significant negative correlation with 
soil clay particle content, cation exchange capacity and pH. Because of the high content of fulvic acid of soil, a positive 
correlation existed between bioavailable soil Pb, Cr and Cd contents and soil organic matter content. Sol bioavailable As 
content was positively correlated with soil pH, but there was no significant correlation between soil bioavailable As and other 
soil physical and chemical properties. The correlation between bioavailability of soil heavy metals and soil basic 
physico-chemical properties was critical for developing appropriate management practices that control heavy metal pollution 
in coastal mudflat regions. The findings of this research provided additional scientific basis for source reduction, ecological 
activity passivation, pollution abatement and risk prevention of soil heavy metals in mudflat coastal areas. 

Keywords: Coastal mudflat; Heavy metal; Spatial distribution; Land use type; Soil physical and chemical properties; Jiangsu 
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Table 1 Descriptive statistics of the available and total contents of soil heavy metals in the coastal region (7-77) 
EE BMA B vi sara PARR aa RAE F 
Heavy Minimum Maximum Mean Coefficient of Natural background Soil environmental quality standard (mgkg ) 
metal (mg:kg )) (mg:kg )) (mg:kg ') variation (96) ind . — Ho ZARATE 
(mgkg ) First class standard Second class standard 
Pb, 9.21 68.3 28.0 61.48 11.40 35 350 
Pb, 0.287 2.30 1.14 54.83 
Cr 17.6 75.9 56.1 21.35 60.11 90 250 
Cr, 0.202 0.578 0.397 20.21 
Cd, 0.027 0 0.103 0.072 7 34.73 0.042 0.20 0.60 
Cd, 0.013 5 0.034 1 0.0214 19.33 
AS, 5.21 9.35 7.17 11.50 7.38 15 25 
AS, 0.041 5 0.123 0.064 3 21.19 


BSeRash finvRRSSSS, Fin‘a’RRAWMS SAH. The subscripts "t" and “a” following symbols of heavy metals indicate the total 


contents and available contents of the heavy metals, respectively. 
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Table 2 Spatial trends of the total and available contents of soil heavy metals in the coastal region 


BEB 48% AM Correlation coefficient ZB RU 

uou 42/8 Longitude (Yi) 46 Latitude (X2) Multiple linear regression 
Pb, -0.774" -0.478" Pb, —4.58E-03X,4.65E-04X5*3 081.647 (r=0.775, F=55.5, P«0.001) 
Cr, 0.445" 0.373" Cr, -1.47E-03.X, -1.17E-03X5-4 632.095 (r=0.468, F=10.35, P<0.001) 
Cd, 0.710" 0.611" Cd, -4.82E-06X,--4.28E-06X5-16.879 (120.752, F=48.08, P<0.001) 
As, -0.181 -0.126 As, —3.92E-05X1-1.26E-05X5-- 64.572 (r=0.183, F=1.28, P=0.285) 
Pb, -0.638"^ -0.371" Pb, — 1.42E-04X,-1.51E-06X5-49.048 (r—0.638, F=25.46, P<0.001) 
Cr, -0.258" 0.072 Cr, —1.28E-05X;-- 1.48E-05X75*49.451 (70.373, F=5.98, P<0.01) 
Cd, —0.047 —0.030 Cd, ——6.63E-08X,—8.556E-09X540.072 (r=0.047, F=0.083, P=0.921) 
Asa 0.063 0.023 As, 73.65E-07X;—1.54E-07X5--0.516 (r=0.066, F=0.160, P=0.853) 


BEBRShhiMvRHSSEs, Paca US A EB ER SUEEZKCE. P«0.05; ** RREA IEKE P<0.01. The subscripts "t" and 
“a” following symbols of heavy metals indicate the total contents and available contents of the heavy metals, respectively. **" means significance at 


0.05 probability level; ***" means significance at 0.01 probability level. 
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Table 3 Normality test result and fitted semivariogram parameters of the total and available contents of soil heavy metals 


ae : (RE ium K-S PTS He Renee We E oss 

meal Skewness Kurtosis P Nugget Sill Nugger/sill (%) Range (m) coefficient (7?) 
Pb, —0.009 -0.705 0.984 10.40 144.60 7.20 3 750 0.861” 
Cr, —0.334 1.195 0.746 45.10 113.80 39.63 2 831 0.698" 
Cd, 0.012 -1.251 0.472 2.80E-05 3.51E-04 7.98 3 920 0.923" 
As, 0.269 0.955 0.865 0.202 0.522 38.70 3316 0.651” 
Pba 0.002 —0.847 0.525 0.010 0.304 3.29 4 470 0.893" 
Cr, —0.640 0.467 0.215 1.95E-03 5.32E-03 36.65 2 675 0.585" 
Cda 1.169 1.310 0.340 3.04E-06 2.34E-05 13.00 3 389 0.839" 
AS, 1.498 4.043 0.178 7.91E-05 1.96E-04 40.36 2 616 0.591" 


BEeBhsh RACCAURÉ$sSS6e, Matan 3US Gi. LIR Pb. Cr. Cd £385 Pb. Cr S US A EECAUA ERG EA UI IBS; + 
J£ As $58 5 Cd, As BASA SHR. The subscripts "t" and “a” following symbols of heavy metals indicate the total contents and available 


contents of the heavy metals, respectively. The residuals with the spatial trend removed were used for the total contents of Pb, Cr, Cd and available 
contents of Pb and Cr; The raw values were used for the total content of As and available contents of Cd and As. 
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Table 4  Cross-validation results of spatial estimation of the total and available contentsof soil heavy metals 


ku 1-37 FRE HARRE EAM HEIR 1375 MAIR 18 X RR 
Soil heavy metal ME RMSE ASE RMSSE Regression coefficient (r) 
Pb, 0.064 7.28 7.10 1.03 0.91" 
Cn 0.81 8.95 8.61 1.06 0.67" 
Cd, 2.64E-04 0.013 0.012 1.03 0.87" 
As, —1.58E-02 0.58 0.55 1.09 0.52" 
Pb, —2.07E-03 0.33 0.29 1.13 0.85" 
Cr, —1.40E-03 0.063 0.066 0.95 0.62" 
Cd, 3.55E-04 3.42E-03 3.71E-03 0.99 0.57" 
ASa 3.85E-04 0.012 0.013 0.95 0.46" 


BSEERESBURAmCAURSESE, Fin‘a’RRAMS EAH The subscripts "t" and “a” following symbols of heavy metals indicate the total 


contents and available contents of the heavy metals, respectively. 


PHR ME(Mean Error), SJARIRZ RMSE 
(Root Mean Square Error) fl PHAR ER ASE(Average 
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Fig.3 Spatial distribution of the total contents of soil heavy metals 
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Fig.4 Spatial distribution of the available contents of soil heavy metals 
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